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We show that pQCD factorization incorporated with pre-haronization energy-loss effect naturally 
leads to flatness of the nuclear modification factor Raa for produced hadrons at high transverse 
momentum pr- We consider two possible scenarios for the pre-hadronization: In scenario 1, the 
produced gluon propagates through dense QCD medium and loses energy. In scenario 2, all gluons 
first decay to quark-antiquark pairs and then each pair loses energy as propagating through the 
medium. We show that the estimates of the energy-loss in these two different models lead to very 
close values and is able to explain the suppression of high-pT hadrons in nucleus-nucleus collisions at 
RHIC. We show that the onset of the flatness of Raa for the produced hadron in central collisions 
at midrapidity is about px ~ 15 and 25 GeV at RHIC and the LHC energies, respectively. We show 
that the smallness {Raa < 0.5 ) and the high-pr flatness of Raa obtained from the kx factorization 
supplemented with the Balitsky-Kovchegov (BK) equation is rather generic and it does not strongly 
depend on the details of the BK solutions. We show that energy-loss effect reduces the nuclear 
modification factor obtained from the kx factorization about 30 -r 50% at moderate px- 



I. INTRODUCTION 



The nuclear modification factor (NMF) Raa is defined as 

d 2 N AA 
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where N co u denotes the number of binary collisions in nucleus-nucleus (AA) collisions. The measurements in Au + Au 
collisions at RHIC showed that the value of Raa is small and constant up to high transverse momentum of produced 
hadrons in central collisions L see Fig.Q] At first sight such a behavior at high-pT contradicts the perturbative QCD 
(pQCD) factorization theorem 7Hll|. Accordingly to the factorization theorem, the inclusive gluonic-jet production 
cross-section with a large transverse momentum px is proportional to ^4 2 < J j iarc ja; G 2 {x,px) which leads to Raa 1 
(see Sec. Ill for the details). On the other hand in the Color Glass Condensate (CGC) approach due to gluon saturation 
and the appearance of a new dimensional scale, saturation momentum Q s 12- 16] , there is a priori no reason to expect 
the pQCD factorization theorem to be valid. Indeed, it has been proven [17h23 | that the factorization theorem can 
be replaced by the kx factorization [23l - [26j for scatterings of dilute-dilute or dilute-dense system of partons (like 
deep inelastic scattering with nuclei or scattering of two virtual -photons). For the case of scatterings of dense-dense 
systems when we have three scales: two saturation momenta for two nuclei and the transverse momentum of produced 
jet, the kx factorization has not been proven yet. However, for understanding of the NMF behavior, it is enough 
to discuss the behavior of the inclusive hadron production at transverse momentum px larger than both saturation 
scales. In this kinematic region, the kx factorization works [l7| . Having this in mind we can assume that the kx 
factorization works in all kinematic regions for scattering of dense-dense systems. 

Based on the kx factorization, the inclusive production for gluon (jet) in AA collisions at midrapidity can be 
calculated from the following equation, 

^^=^krkS dHd2B ^o(^-b\), ( 2) 

where px and y are the transverse momentum and rapidity of the produced gluon, with notations z = ln(r 2 Q 2 ), 
x± = Pt/Qs and Nq — 2N — N 2 . We defined Cf = (N 2 — 1)/2N C where N c denotes the number of colors and a s 
is the strong coupling. The forward dipole-nucleus amplitude N can be obtained via solving the Balitsky-Kovchegov 
equation (BK) 15]. It is seen that Eq. (J2|) has a^-scaling behavior in the kinematic region that the forward dipole 
amplitude N has the geometric-scaling property. 

Using Eq. |2j) we can re-write the NMF in the following form: 
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FIG. 1: The NMF for hadrons and direct photon versus transverse momentum pr in central Au+Au collisions at RHIC. The 
plot is taken from Ref. [f| from the PHENIX collaboration. 



where function T is defined using Eq. © . The area of interaction (denoted by Sa and S p for AA an d pp collisions) 
is well-defined for nucleus-nucleus scatterings while is rather uncertain for proton-proton interactions |27l [28j . In the 
above for simplicity we assumed that the area of interaction in both AA and and pp collisions can be factorized. 

It turns out that the NMF calculated from Eq. (J3]) is small Raa < 0.5 (see Sec IV). Eq. ^ has the geometric scaling 
behavior which is valid in the saturation region [29j and can be extended to a wider region outside the saturation 
domain [3(|. The geometric scaling behavior is violated for z 2 > 2 In [Q 2 {Y) / 'Q s (0)) where Y — ln(l/x) [3(j. Notice 
that wc have 2 In (Qf (Y)/Q a (0)) instead of Ehi(Q 2 (Y)/Q s (0)) since d 2 a/d 2 p t oc oc iV 4 when N is not small. 
Therefore, for z ^> y/2 In (Q 2 {Y) / Q s ({))) we expect that the scaling behavior will be broken and the inclusive cross 
section for jet production for both nucleus- nucleus and proton-proton collisions will be proportional to a s {p 2 ^Jp T 
and consequently we have Raa 1- Using the above equation and the saturation scale estimated in Ref. [3l| at 
RHIC energies, one may expect that only for px < 3 -j- 4 Q s , the NMF to be different from unity. But a slight glance 
at Fig. Q] shows that this is not the case for inclusive hadrons production for a wide range of px while it apparently 
holds for the production of direct photon [32l434l | . Therefore, the small value of the NMF at high px may stem from 
larg e distance interactions in medium and non-perturbative QCD should be invoked to calculate such interactions 
(35j- The data on Raa for J/* production supports this conclusion j37l - l40l |. 

In this paper, we show that the CGC predicts considerable suppression of the NMF at low-pr both at RHIC and 
the LHC energies. We will show that based on the CGC prescription at low-py, the NMF slowly increases and then 
flattens at moderate pr- Based on pQCD, we will then estimate, the onset of the flatness of the NMF at large 
pr- The CGC approach is based on classical gluon fields and evolution in leading log(l/x) approximation. In such 
approximations the energy-loss in processes of gluon emissions is neglected. The systematic approach within pQCD 
framework [3o1 . l4l| - |46| shows that such emissions alone is not able to describe the experimental values for the NMF and 
variety of non-perturbative approaches have been developed (46l - [54| . In this paper, we introduce a simple approach 
for calculating the energy-loss effect due to pre-hadronization which is able to describe RHIC data at high-px- At 
Icm-pT, energy-loss is less important and we show that the CGC prescription Eq. Q gives rather good description of 
RHIC data at small Bjorken-x. 

Our approach is based on two main ingredients: (i) the kx factorization [l7l - r22j | which takes into account the 
rescatterings of the produced dipole (or gluon) in the final state; and (ii) the energy loss due to the hadronization of 
the produced dipole (or gluon). 

In the next section, we calculate the energy-loss effect within two different pre-hadronization pictures. In Sec. Ill, 
we discuss the large px behavior of the NMF for hadron production at RHIC and the LHC. Sec. IV is devoted to 
present our numerical results and comparison with the experimental data. As a conclusion, in Sec. IV we highlight 
our main results. 
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II. PRE-HADRONIZATION MODELS 

Hadronization, unfortunately, could be treated mostly phenomenologically due to our lack of understanding of 
non-perturbative QCD. There have been several approaches to describe hadronization processes [46 54]. 

Here, we propose a simple picture for hadronization hoping that the phenomenological uncertainties can be reduced 
to few parameters which can be then determined from experiment. In our approach, a high energy heavy ion collision 
can be simplified to three successive stages: 1) parton production 2) pre-hadronization and 3) fragmentation of the 
produced parton into hadron. At the first stage, we have partons scattering and recombination effects which can be 
described within the CGC approach. The CGC gives the description of the initial wavefunction of projectiles and 
the multiple-parton production within a time scale of the order of 1/Q S after collisions. In the second stage, the 
produced dipoles propagate, emit gluons and interact with the produced QCD medium and lose energy. As we have 
already mentioned, the emission of gluons and interactions with other dipoles have been taken into account within 
the kr factorization if the sizes of dipoles are small enough to trust the perturbative QCD approach. In the CGC 
approach most of produced dipoles have sizes of the order of 1/Q S <C where // is the soft scale. The energy loss 
for such dipoles can be taken into account by the BK equation with a kernel that includes the energy conservation. 
The estimates of Ref. [65| shows that these corrections cannot explain the experimental data leading only to 20-^30% 
suppression of the NMF. We believe that the main source of the energy loss at this stage of process is originated from 
hadronization of the slow dipoles that can be created during interaction of the fast dipole with the slow remnants 
of the nuclear target or with the fast remnants of the projectile nucleus. In the lab. frame one of the nucleus is at 
rest and the produced dipole can interact with a nucleon in the target or projectile elastically or inelastically. Elastic 
interaction does not break the nucleon and does no lead to any energy loss. However, in inelastic interactions nucleons 
are broken, and due to hadronization, the remnants of the nucleon carry (1 — z^) fraction of energy leading to z^ 
energy loss of the fast dipole. The produced hadrons do not interact since the density of the nucleons at fixed impact 
parameters are small in the nuclei (maximum 2 nucleons for gold ). Finally, the sizes of produced dipoles during 
interaction with nuclei reaches the soft scale (1/m) and due to hadronization the energy of the fastest hadron becomes 
Zh times smaller. Notice, that we assume for simplicity that both hadronization stages have the same Zh 

In our description of hadronization, we neglect the fragmentation process and instead for simplicity and clarity 
resort to the Local Parton-Hadron Duality (LPHD) principle [13, [H, [3l|, H|| , namely we assume that the final-state 
hadronization is a soft process and cannot change the direction of the emitted radiation 1 . Then, the key question 
is to obtain the fraction of energy of the produced partonic system (in the pre-hadronization stage) which is carried 
away by hadron, the so-called ZA,h parameter. We have recently shown [27| that the data for the inclusive hadron 
production in pp collisions in a wide range of energies can be described if we assume that the energy loss in the process 
of fragmentation of the produced gluon into hadron reaches 2/, w s, In the following section we calculate the value 
of za in the presence of a dense QCD medium within two different pre-hadronization pictures and we will show that 
both schemes give rather similar results. 

A. Quark-antiquark pre-hadronization model 

As a first model, in spirit of all available hadronization models, we assume that hadronization passes a stage in 
which all gluons first go to quark-antiquark pairs (pre-hadronization) (see Fig. [2]). Each pair then decays into hadron 
in the same way for hadron-hadron, hadron-nucleus and nucleus-nucleus interactions. In such a decay, the fastest 
hadron carries Zh fraction of energy of the pair. 

In nucleus-nucleus interactions, the process of quark-antiquark pair creation is depicted in Fig. [3J In heavy ion 
collisions, the probability for a gluon 2 with an average transverse momentum Q s \ to decay into quark-antiquark (qq) 



1 The same idea was used in the KLN [3ll and the LR [2H approaches which describes the rapidity distribution of heavy-ion collisions 
data in a wide range of energies. This scheme also describes experimental data from e+e - annihilation into hadrons 15511 and inclusive 
hadron productions in pp collisions including the recent LHC data |27H . 

2 It is well known that at high energies the most produced partons are gluons. However, at RHIC energies the contribution of quarks, 
perhaps, are not negligible. However, in our approach similar to Ref. [36] , we assume that gluons give the most important contribution 
at the highest RHIC energy as well. The main reasons behind this assumption are: the main contribution for the NMF factor measured 
at RHIC, stems from the momentum of the jet pj- = 8 — 20GeV for z = 1/2. The typical values of x = 2px/W/z turns out to be 
smaller than 0.1 for W = 200 GeV while the quarks become essential at x > 0.2. Moreover, our main goal here is not to describe 
the experimental data but to answer the question whether it is possible to have the value of the NMF much smaller that 1 within the 
framework of the QCD factorization theorem. 
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pair at impact-parameter b with a relative transverse momentum denoted by k = k (fc 2 = k 2 ) has the following form 3 
P q A(v^,k) = ^0- J d 2 r J d'r'e-^-^^^KtirQ^Ktir'Q^+KoirQ^Koir'Qsi)} 

~ J^y (l " e-i^"*) (l - e-i(«3'*) } . (4) 

The factor in the first curly bracket is the distribution function ^(G — > qq(r)^*(G — > qq{r') that describes the 
production of qq dipoles by gluon with size r (r') in the amplitude and complex conjugated ampplitude, respectively, 
before the passage through the nuclei (see Fig. [3]). This gluon has typical transverse momentum Q s \. The formula is 
written in the laboratory frame where the nucleus Ai is at rest. Eq. (0| is derived as a sum of all possible inelastic 
interactions. For example, for one inelastic interaction with nucleus A-i the contribution to Eq. (Q} has the following 
form, 

dzie~? a ( Y ' r)pzl pa in e^ a{ - Y ' r)p{2RA ^ Zl) e~^ Y ' r) p2Ra i e ^^ a{Yy) p(2Ra ^ Zi) e^^' j(Y ' r " ) pzi , (5) 

where p is the nucleon density and R Ai is the radius of the nucleus Ai. Eq. ([5]) indicates that a dipole with a transverse 
size r interacts inelastically with the nucleon located at point X\ in nucleus A 2 and exponents provide that no inelastic 
interactions occur with other nucleons in both nuclei in the amplitude. In the complex conjugate amplitude the same 
happens with a dipole with a transverse size r'. Note that for simplicity we assumed that the nuclear profile is 
cylindrical T(b) « 2R A with b 2 < R\. 

The dipole-nucleon cross-section can be written as 

47r 2 a s f d 2 l ,„ ir.i\ -tvi\ ,( V i2\ rl«i 2TT 2 a s 2 



a (Y,r) = ^ J (1 _ e *trL) (i _ e -*vL) * (y, ^ ^ r 2 J ^ (y, , ( 6 ) 

where is the unintegrated gluon density, y is the rapidity of the dipole. In the right side of the above equation 
we have kept only the dominant contribution relevant at large transverse momentum of the produced dipole in the 
DGLAP approximation. Using Eq. ^ , one can simplify Eq. (O by the following replacement [H, [H[ : 

a(Y,r) pT(b)^a(Y,r 2 ) p2R A -> r 2 Q 2 sl /8, (7) 

where Q S i is the gluon saturation momentum in the nuclei i and can also depend on the impact-parameter b. In the 
same fashion, for inelastic dipole-nucleon cross-section we have the following expression 

^(Y,r,r') = 2^ j ^ (l - 4*1) (i _ e -^i) , ( y ,») ^ ^ 2 , V j dl^{Y^) , 

a m (y,r,r') p2i? A ^2r-/Q 2 l /8. (8) 

The inelastic cross-section can be generally written in a form of convolution of two amplitudes <Ji n — ^Z n A(r;{ki\)x 
A* (r'; {hi}) with sum over produced n-partons with momenta ki. Eq. ([8]) describes the inelastic cross-section for 
the process in which the dipole of size r in the amplitude interacts with nucleon while in the complex-conjugated 
amplitude, the dipole with size r' interacts with the same nucleon as in the amplitude. 

Note that Eq. Q integrated over the impact-parameter b gives the inclusive cross-section for gg-pair with a relative 
transverse momentum k. It incorporates the first stage of the process, namely the parton production (shown as a 
dark blob in Fig. [5]). Here, we assumed that the produced quark-antiquark pair is moving in the classical gluon field, 
neglecting the evolution in the rapidity interval y. We believe that this assumption will not be essential in estimating 
the energy loss. 

In calculating the average energy loss, we assume that in each inelastic interaction, the dipole energy decreases 
Zh times due to the pre-hadronization process. For elastic interaction we also introduce factor z\ lfi i which describes 
the fraction of energy that is carried by the fastest hadron. For simplicity we assume that Zh — Zh, e i- For elastic 



3 The derivation is based on the approaches suggested in Refs. 17, 33, 59-61]. 
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FIG. 2: Inclusive production of hadrons in nucleus-nucleus collision through the stage of pre-hadronization. 




FIG. 3: Inclusive production of quark-antiquark pair (dipole) in nucleus-nucleus collisions. 



scattering the energy loss emerges only on the last stage, namely quark-antiquark (or gluon) hadron transition. An 
additional source of the energy loss in nucleus collisions stems from the inelastic interaction of gg-pairs (or gluons). 
In this process, soft gluons (slower than propagating qq-paiv or gluon) are produced. Then, these gluons produce 
hadrons due to hadronization and these hadrons carry (I — Zh) energy of colliding gluon. Due to this process, the 
propagating gluon loses Zh fraction of energy per each inelastic interaction. Therefore, in order to calculate the energy 
loss in a nucleus we need to find the probability to have the inelastic collisions with i nuclcons in a nucleus (Pi) and 
multiply it by the energy loss in each such collisions. In other words, 

<z A > = z h P + ^2 zi p h (9) 
i=l 

where we define the probability Pi = (Ji,i n /&in with Oi^ n being the contribution of inelastic i collisions to the total 
inelastic cross section <7j n . The latter is given by Eqs. (|4l5p . Therefore, in order to calculate the resulting average 
energy- loss of the dipole propagating through the medium, in Eqs. (I4I5|) . one needs to find a^m and multiply every 
inelastic collision by the factor z^. From Eqs. (I6I8[) one can see that each factor r ■ f_ is associated with the inelastic 
interaction. Therefore, this multiplication can be carried out by the following replacement in Eq. 

(r - L_f = r 2 + r' 2 - 2r ■ r_ ->■ r 2 + r' 2 - 2 z h r ■ r_. (10) 

The fact that factor r • rf_ is associated with the inelastic cross-section in Eq. ([8]) for Pj^a (d> ^ ^) can be understood 
on general grounds. For example, the inelastic dipole cross-section at two-gluon exchange level is Gi n oc rr' . Notice 
that all other terms in expansion of P\ 9 A (y, b, k) are of the order of r 2 and r' 2 . 
The average energy loss can be then obtained from the following expression 

(z h ) A = z A = ((z h - 1) Pf A (y, b, k- z h = 0) + Pf A (y, b, k- z h ) ) / Pf A (y, b, k; z h = l). (11) 
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In the above equation, the term P'aa (y, b, k; Zh = 0) describes the process of diffraction dissociation in which dipoles 
scatter only elastically and the term P^a (v? b, k; Zh) — P'aa iVi ^> ^> Zh ~ 0) corresponds to the contribution of all 
inelastic processes with additional factor Zh for each rescattering. 

Therefore, the first term in Eq. only takes into account the processes of hadronization of the quark-antiquark 
with rapidity y that passes through the medium without an inelastic interaction. The second term in Eq. (1111) is 
responsible for hadrons production in the entire kinematic region of rapidity y. Note that Eq. (jlll) is based on the 
assumption that if no energy-loss occurs in a single interaction Zh = 1 then no energy-loss will follow in multiple- 
interactions namely Za = 1, and if the parton loses all its energy in a single interaction — then naturally in 
multiple- interactions we have also za — 0. The NMF Eq. ([3]) incorporating the pre- hadronization energy- loss effect, 
can be then rewritten in the following form, 

1 Si T[pt / {za Qs.a) J 

A S P T\Pt J ( z h Qs,n) J 

In Eq. (|T2j) we use the fact that the transverse momentum of the emitted gluon in pp or AA collisions is equal to [13] , 

p\ (gluon) = p\ (hadron) /z h A + fcf (intrinsic) , (13) 

where kx (intrinsic) is the average intrinsic transverse momentum of the gluonic jet and since its value is small, it can 
be neglected. 



B. Gluon pre-hadronization model 

In the second model, we consider a case that the produced gluon decays into hadrons after propagating through 
the nucleus (see Fig. |4]). 

Using the approaches suggested in Refs. [lj], HiJ [6(| and summing all gluon exchange diagrams shown in Fig. [S] 
and neglecting again the small- a; evolution in the rapidity interval between the produced gluon and the target, the 
probability of gluon production in AA collisions with a transverse momentum fc can be written in the following form, 



P£tM,A0 - ^f 1 J d 2 r J dVe-^-^ \k 2 (kr) K 2 {kr') J (kr) J (kr') j 

- J^y (l - e"*<«^) (l - e-*(^"*) } . (14) 

In the derivation of the above equation, we assume that the gluon with transverse momentum k is emitted by the 
dipole with transverse size smaller than 1/k. Indeed, all rescattering of the gluons emitted by dipoles with larger 
sizes are already included in Eq. (j5J| 17]. In this case the emission of gluon from interaction with a nucleon can be 
considered as a decay of the nucleon to two colorless dipoles. This can be better conceived at large N c limit where 
an adjoint (gluon) dipole can be decomposed into two fundamental (quark) dipoles. These two dipoles penetrate 
through the medium with the same cross section. The wave function of such a gluon has been calculated in Ref. 62] ■ 
Therefore, this approach differs from the quark-antiquark pre-hadronization model mainly by the value of the dipole 
cross section: oq (Y,r 2 ) » 2<Jdi P oi e - In other words, the main differences between Eq. (Ti"4| and Eq. (|4]) is that in 
Eq. (TTJ} the value of the saturation momentum is two times larger than in Eq. 

Notice than in general, we are dealing here with rescatterings of two dipoles with different sizes. For example, in 
the dipole-nucleus scattering amplitude at the first stage of the process, the incoming dipole with size X12 — |r*i — f%\ 
emits a gluon with a coordinate r and this process at N c ^> 1 can be considered as a decay of the dipole X12 to two 
dipoles X13 — |fi — r\ and X23 = \r2 — r\. These two dipoles take part in interactions with a nucleus. However, in 
the inclusive cross-section we integrate over the momenta that carry quark (antiquark) with coordinates r\ and r 2 - 
Therefore, this cross-section does not depend on these coordinates and both can be put equal to zero. That is the 
reason why the inclusive cross-section here can be viewed as a propagation of dipole with size r in the amplitude and 
with size r' in the complex conjugated one. 

In order to calculate the average energy loss of the produced gluon at the pre-hadronization stage (zh)A, we again 
follow the steps in Eqs. (|10lll[) by replacing P'aa ~^ ^aa- The corresponding NMF can be then obtained by Eq. (|12l) . 
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FIG. 4: Inclusive production of hadrons in nucleus-nucleus collision due to gluon decay. 




FIG. 5: The passage of the produced gluon in nucleus-nucleus collision. 



It is worthwhile mentioning that in Eq. (|12l) we take into account both the first stage, namely gq-pair (or gluon) 
production (the dark blobs in Figs. [2JH]), and also the energy-loss effect due to the passage of the produced gg-pair (or 
gluon) through the medium (see Figs. [31 [SJ. In a sense, we unfolded the kr factorization and made the corresponding 
corrections due to the pre-hadronization energy-loss effect since the hadronization leads to its violation. 



III. LARGE p T BEHAVIOR OF NMF 



At very large pt, the pQCD factorization theorem can be used. The production of gluon-jet with transverse 

momentum px (sec Fig. [5]) in nucleus-nucleus collisions can be simply written as 

dy?dPpr ^ y=0 = A * d^dPpr = ah ^ dXlGA = 2pT ^' pT ) X2 ° A ( X2 = Wv 7 ^, ~Pr) , (15) 

= A 2 ^± Xl G p ( Xl = 2 PT /^p T ) x 2 G p (x 2 = 2 PT /VS, -pt) -+ A 2 ^M (p^/Qlf , 
Pt Pt 

where xGp (or xGa) is the gluon structure function for the proton (or nucleus) and Qo is & separation scale so 
as for pt > Qo one can use perturbative QCD. In the above expression, cr narc [ is the perturbative partonic cross 
section computable up to a given order in a s , and 7 denotes the anomalous dimension. It is well-known that in the 
leading-order for x — > 1 we have 7 — > a s . Therefore, at very large pt 3> Q s but pt -C the NMF obtained via 
Eqs. (|1I13I15|) is equal to 



a 2 {pT/z h )(pT/z h f as(PT/Zh) 

The value of the NMF can be then immediately estimated at large pt via Eq. (fTB| by only knowing the value of Za, 
Zh and also Aqcd which appears in the running strong-coupling. For the running coupling a s , we employ the scheme 
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FIG. 6: The gluon production accordingly to the factorization theorem at large pr- 



used in Ref. [27[ at the leading-order. We postpone the numerical discussion to the next section, and in the rest of 
this section, we address another important question, namely at which px the expression given in Eq. (|16l) can be used. 
Let us for sake of simplicity work in double log approximation with the anomalous dimension 7 = a s /w where uj is 
Mellin image of ln(l/x) and ct s = a s N c /Tr . In such an approximation, the nuclear gluon structure function up to a 
numerical constant JC can be written as 63] 



xG A (x = 2p T /y/s,p T ) 



Pt 



/C exp J4a s ln(l/z) In (p|/Q2) - In (p 2 T /Q 2 ) + l A , 



= /Cexp [JDntpl/QD - MPt/Qo) + U ), 



(17) 



where L = 4a s ln(l/a;) = In (Q 2 (A; x) /Q 2 S (A; x = xo)) and Ia = (1/3) In A. The equation for the saturation scale 
has the form 1631. 



L l s = (l s — Ia) 2 



with a solution, 



l s (L, l A ) = In (Q 2 (A; x) /Q 2 ) = ^+l A J(^) + L I x . 



(18) 



(19) 



One can see that Ga given by Eq. (fT7|) at Q 2 = Qq is proportional to A 1 / 3 . Introducing a new variable Z = 
In (pt/Q 2 (A; x)) one can rewrite Eq. (ITTl) in the following form 



xG A (x = 2p T /y/s,p T ) 
Pt 



K exp(H(Z,L,l A )) , 



with a notation, 



H(Z,L,l A ) = VL {Z + l s (ZJa)) - Z - l s (L,l A ) + l A . 
Expanding H (z, L, Ia) at small z, we have 

H(Z,L,l A ) = -h 1 (L 1 l A )Z + h 2 {L,l A )Z 2 + 0(Z 3 ), 

with 

L 



hi {L,l A ) 



h 2 (L,l A ) 



1 



2\l ^ + Ll A 



2\/2 [L [L + 21 A + \/L 2 + ALlj 



■81'. 



3/2 • 



(20) 

(21) 
(22) 

(23) 

(24) 



The function hi and h% versus L for the gold are shown in Fig. [3 It is seen that at large values of L, we have hi 1/2 
and h?, — > 1/ 8L in accordance with finding in Ref. [30I ] . Using the KLN saturation model parametrization, the value 
of L is about L«l and 1.7 for RHIC and the LHC energies, respectively [31|. The condition [30J 



h 2 Z 2 > 1, 



(25) 
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FIG. 7: Function hi and h 2 denned in Eqs. (|23l24j) versus L. 



characterizes that Z is large and the scattering amplitude is far away from the saturation domain, out of the region 
where we have the geometric scaling behavior at large values of L. Indeed, the first term in Eq. (|22[) leads to the 
geometric scaling behavior at large values of L while the second term violates this behavior even at large L. Therefore, 
the condition given in Eq. (f^5|) ensures that the gluon structure function defined in Eq. (f2T))) gives the perturbative 
form used in Eq. (IT51) . 

Eq. ([25]) leads to pr > 40 GeV for the RHIC energies and to pr > 60 GeV for the LHC energies for the produced 
gluon. In terms of the transverse momenta of the produced hadron we expect that for p? > 40 za 15 GeV at the 
RHIC and for p T > 60 z A « 25 GeV at the LHC, the NMF defined in Eq. (16]) to be reliable. Notice that if the 
pre-hadronization stage in A A and pp collisions was the same, namely za = Zh, then the NMF defined in Eq. (|16[) was 
identically equal to one Raa = 1- Therefore, the main source of suppression of the NMF for hadrons at high-px is 
due to further energy loss of the produced gluons in AA collisions compared to pp collisions at the pre-hadronization 
stage, namely za < z p . Notice that on average za and z p vary slowly with px for the kinematic region of our interest 
here [2g|. Note also that the pT-dependence in Eq. (TTB1) mainly enters through the running of the strong-coupling. 
This justify the flatness of the NMF Raa over large range of pt- In the next section we show that the numerical value 
obtained for za from the pre-hadronization mechanisms introduced in section II can indeed describe the observed 
suppression of hadrons at large-pT at RHIC. For smaller values of pt, Eq. ^ may be used with the dipole amplitude 
N being in the region where we still have the geometric scaling behavior for N = N(z) (30j . 



IV. NUMERICAL RESULTS AND COMPARISON WITH EXPERIMENTAL DATA 

In order to estimate the value of {zhjA, we use the KLN saturation parametrization [3lJ with Q s — 1.41 GeV at 
\fs = 200 GeV for gold and we take k = 2 GeV as a typical scale for soft interaction. The chosen soft interaction 
scale corresponds to the value of the soft Pomeron slope k « 1/ J dp with a' P — 0.25 GeV -2 [6(1 HH, |64j. It turns 
out that both pre-hadronization models give very close values for (zh) a/ Zh ■ 0.766 and 0.744 for the quark-antiquark 
and the gluon pre-hadronization model respectively, at the RHIC energy y/s — 200 GeV. The value of Zh ~ 1/2 was 
obtained from a fit to the ex per imental data in pp collisions for the average transverse momentum of charged hadrons 
for a wide range of energies [271 ] . 

Again using the KLN parametrization for the nuclear saturation scale, we estimated the value of (zh)A/zh at the 
LHC energies, assuming Zh ~ 1/2. At higher energies or larger value of the saturation scale Q s , one can approximate 
the integrand of Eqs. (|4IL4[) to its maximum value at r = rl w 1/Q S . This leads to a similar value for za in both 
pre-hadronization models at higher energies, since the over-all factors in Eqs. f|4IL4[) will be canceled out via Eq. (|lll) . 
Indeed, full numerical solution is in accordance with this observation and we found that both pre-hadronization models 
result to a very similar value for za at the LHC energies, namely (z^a/zh ~ 0.82, 0.795 at ^/s = 5.5, 2.75 TeV, 
respectively. 

Our results for the NMF for — 10% centrality in Au + Au collisions using the k t factorization Eq. (fl~2"|) supplemented 
with solutions of the BK and the modified BK equation that preserves the energy conservation [65J is shown in Fig. HI 
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FIG. 8: The NMF for hadrons obtained by Eq. (|12|) for — 10% centrality Au + Au collisions at various energies. Function 
T(x±) is calculated using the solution of the BK and the modified BK equations. The solid and dashed lines correspond to 
the quark-antiquark and gluon pre-hadronization models, respectively. The experimental data is taken from Ref. PJ. 



In the modified BK equation (denoted by ModBK in Fig. the next-to- leading order corrections to the BFKL 
kernel are taken into account which leads to conservation of energy in the framework of non- linear equation 66] . We 
recall that here we resort to the LPHD principle for the final state hadronization. But the energy-loss effect in the 
pre-hadronization both in pp and AA collisions are effectively incorporated in the results shown in Fig. [5] We believe 
that the difference between two pre-hadronization models (introduced in Sec. II) in AA collisions gives reasonable 
estimates of the errors which are unavoidable in the situation that we do not know the theory of pre-hadronization 
in a dense medium. As we already pointed out, the differences between these two pre-hadronization models at LHC 
is negligible. The energy-loss effect reduces the NMF defined in Eq. (|12p about 30 -j- 50%. Notice that this reduction 
does not depend on a given BK solution and incorporates the missing energy-loss effect into the kx factorization. 

We should stress that the numerical estimates in Fig. [8] was obtained from the BK equation with a simplified kernel 
[63l [65[ . Such estimates is less reliable for describing the experimental data, even though, interesting enough, such 
a rough approximation apparently is in agreement with RHIC data at high-p^ within the errors. Notice that in fact 
the validity of the kx factorization Eq. (fT2j) beyond the extended geometric-scaling region at high-pr is questionable, 
see also Sec. I. An important message here is that the smallness and the high-pT flatness of the NMF obtained from 
Eq. (|T2"j) and the BK equation is rather generic and it does not strongly depend on the details of the BK solutions in 
our interested kinematic region here, see also recent interesting paper by Albacete and Marquet [67| . It is also seen 
from Fig. [8] that as one may expect, the onset of the flatness of the NMF at the LHC energy occurs at higher px- 
Fig. [5] indicates that the energy-conservation constrain on the BK equation which is important at higher px tends 
to enhance the NMF. This is due to the fact that the ModBK solution correctly incorporates the behavior of the 
anomalous dimension at large Bjorken-x by matching to the DGLAP value. 

In order to make more reliable prediction for the NMF for the neutral pion tt°, we next re-calculate the NMF 
from Eq. p^|) by using the impact-parameter dependent CGC saturation model in pp (2?J and AA collision [28| . It 
has been already shown that such a scheme describes HERA data at small- x [68[ and the hadron multiplicities in 
nucleus- nucleus collisions at RHIC (28j j and provided correct predictions 27 1 for the inclusive hadron production at 
the LHC. The importance of the inclusion of the impact-parameter dependence in the saturation models and the BK 
equation has been addressed in Refs. [13, [H, [6{|. Our results at RHIC and the LHC energies in — 10% centrality 
Au + Au (at RHIC) and Pb + Pb (at the LHC) collisions, is shown in Fig. |H1 We recall that Eq. (TT2"j) has xj_-scaling 
property which is valid only at low transverse momentum roughly for p < 3 ~ 4Qs- In the saturation model [68] used 
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FIG. 9: The NMF for hadrons obtained by Eq. JT2} at low p T and Eq. {T6J at high p T for - 
at various energies. Function T(x±) is calculated within the saturation (CGC) approach [27l. |2S 
obtained in Sec. III. The experimental data is taken from Ref. 



10% centrality AA collisions 
. We also show pQCD results 



in Fig. [SJ the saturation scale varies very slowly with energy and for px = 1 GeV at midrapidity for central collisions 
we have Q s = 0.76 GeV for proton at y/s — 5.5 TeV. Therefore, the upper limit of xj_-scaling behavior in this model 
is about px ~ 3 GeV for pp collisions (the reference for the NMF) at y/s = 5.5 TeV at midrapidity. We therefore, only 
show the CGC results coming from Eq. (|12p for px < 4 GeV. Here, we only concentrate on the CGC results at rather 
low-pr where we can also give a good description of spectra in pp collisions [27]. At the RHIC energy yfs — 200 GeV 
at mid-rapidity, for px > 2 GeV, we have x > 0.01 which is beyond the validity of the CGC prescription. Note also 
that the dipole model used here was obtained from a fit to HERA data for virtuality Q 2 e [0.25,45] GeV 2 . Therefore, 
this model is valid for px € [0.7, 6.7] GeV. Again, we employ the LPHD principle to connect pre-hadronization to 
final hadronization stage. This leads to two different over-all factors for the spectra in pp and AA collisions which are 
fixed with the experimental data at lower energies [13, [HI . The appearance of two different over-all factors for the 
normalization is partly due to the fact that pre-hadronization leads to different effective masses for the mini-jet in AA 
and pp collisions [27l[28j. It is seen from Fig.|5]that our results without any free-parameter to adjust, is in accordance 
with the PHENIX data at low-pT, in contrast to the impact-parameter independent BK solution shown in Fig. [9] with 
a zero mini-jet mass. At \ow-px, the energy-loss in medium is less important and we have taken za = Zh = 1/2 for 
Pt < 4 GeV. Notice that the actual value of Raa at low p T calculated via the k T factorization may be subjected to 
about 15 -i- 25% errors due to the uncertainties in the hadron multiplicity obtained within the same framework in AA 
collisions. Moreover, here the contribution of quarks was neglected which is perhaps less justifiable at midrapidity at 
low-energy. 

At high-pT, we rely on the improved pQCD factorization result given in Eq. (|16|) (also shown in Fig. |9|). The only 
external parameter in the NMF defined in Eq. (fT5]) is the parameter za which incorporates the energy-loss effect at 
the pre-hadronization stage. The value of the parameter za was already calculated in Sec. II (see also above) at 
RHIC and the LHC energies. As we already argued Eq. fp~6|) is valid at high transverse momentum of the produced 
hadrons for px > 40 z a ~ 15 GeV at RHIC and for px > 60 za ~ 25 GeV at the LHC energies. Our predictions for 
the NMF suppression and the onset of the flatness of the NMF at high-p^ both agree with the PHENIX data. As 
we already pointed out Eq. (|T^)) naturally leads to the flatness of the NMF at high px since px dependence mainly 
enters via the running strong-coupling and varies very slowly with px- In Fig. we also show our predictions both 
at low and high px at the LHC energies yfs = 2.75 and 5.5 TeV. The uncertainties in our formulation is less than 
10% which come mainly from the normalization in Eq. (fT2]) . The band in Fig. [9] indicates about 2% theoretical error 
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which also includes the discrepancies between two pre-hadronization schemes. The predictions of other approaches at 
the LHC can be found in Ref. [z^. 



V. CONCLUSIONS 



In this paper first we developed a simple picture for calculating the energy- loss effect. In our approach, the energy- 
loss effect stems from the process of the pre-hadronization. We showed that the estimates of the energy-loss in two 
different models of the pre-hadronization lead to very close values and is able to explain the measured Raa at high 
transverse momentum of produced hadrons at RHIC. We showed that the small value of the nuclear modification 
factor at high-pr is mainly due to the energy-loss effect in the pre-hadronization stage. 

We also investigated the NMF obtained from the kr factorization supplemented with solutions of the BK and the 
modified BK equation that preserves the energy conservation. We showed that the smallness (Raa < 0.5 ) and the 
high-pr flatness of the NMF obtained from Eq. (fl"2j) and the BK equation is rather generic and it does not strongly 
depend on the details of the BK solutions. We showed that the modified BK solution which includes the energy 
conservation tends to slightly enhance the NMF at high-p^- This is due to the fact that the modified BK equation 
properly includes the anomalous dimension at x « 1 which coincides with the anomalous dimension of the DGLAP 
equation while the BK equation only reproduces its double log limit. This indicates that one should be cautious as 
extrapolating the hp factorization results to higher pp is less reliable if the running anomalous dimension obtained 
from the BK solution does not match to the DGLAP value at high-pr- 

We showed that at high-p^, the pQCD factorization incorporated with the pre-haronization energy- loss effect 
Eq. (fT6|) results naturally to flatness of the NMF. We obtained the onset of the applicability of Eq. (fl~6|) to be for 
Pt > 40 GeV at RHIC and for pr > 60 GeV at the LHC for the produced gluon. Notice that in this kinematic region 
the NMF for the produced jet should approach to one R%a 1- However, the NMF of the produced hadron is very 
small R\ A < 0.5 due to the pre-hadronization effect, namely za/z p < 1. In terms of the transverse momenta of the 
produced hadron we expect that this to be for pr > 15 GeV at RHIC and for pr > 25 GeV at the LHC. As one can 
see from Figs. [H [9] the hadronization model is responsible for 30 ~ 50% of the NMF value while 50 -j- 70% of this value 
stems from the BK equation and kp factorization. 

One of important check of our hadronization model will be the estimates of the NMF in deuteron-nucleus collisions 
which will be published elsewhere. Here we would like only to mention that in our approach, estimates for the NMF 
in proton-gold collision at the RHIC energy \/s = 200 GeV at midrapidity gives only upto 10% contribution due to 
the hadronization leading to the value of the NMF about 0.95 ± 0.05. 

Finally, it is worth mentioning that the main assumption that in the process of hadronization the fastest particle 
carries the final part of the initial energy is a typical feature of all hadronic reactions as well as the production of 
hadrons in DIS. 
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